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ABSTRACT 

We present the results of the analysis of the broad-band spectrum of Cygnus X-l 
from 3.0 to 200 keV, using data from a lOksec observation by the Rossi X-ray Timing 
Explorer. The spectrum can be well described phenomenologically by an exponentially 
cut-off power law with a photon index T = 1.451q 02 ( a va l ue considerably harder 
than typically found), e-folding energy E f — 162lgkeV, plus a deviation from a power 
law that formally can be modeled as a thermal blackbody with temperature kT BB = 
1.2^0 jkeV. Although the 3-30 keV portion of the spectrum can be fit with a reflected 
power law with T = 1.81 ± 0.01 and covering fraction / = 0.35 ± 0.02, the quality 
of the fit is significantly reduced when the HEXTE data in the 30-100 keV range 
is included, as there is no observed hardening in the power law within this energy 
range. As a physical description of this system, we apply the accretion disc corona 
models of Dove, Wilms & Begelman (1997a) — where the temperature of the corona is 
determined self-consistently. A spherical corona with a total optical depth r = 1.6 ±0.1 
and an average temperature kT c = 87 ± 5keV, surrounded by an exterior cold disc, 
does provide a good description of the data (x? e d = 1*55). These models deviate from 
the data by up to 7% in the 5 — lOkeV range, and we discuss possible reasons for these 
discrepancies. However, considering how successfully the spherical corona reproduces 
the 10 — 200 keV data, such “photon-starved” coronal geometries seem very promising 
for explaining the accretion processes of Cygnus X-l. 

Key words: radiation mechanisms: non-thermal - radiative transfer - X-rays: bina- 
ries - accretion, accretion discs 


1 INTRODUCTION 

Cygnus X-l is one of the most firmly established persistent 
galactic black hole candidates (BHCs). Since it is also one 
of the brightest BHCs, the study of its X-ray spectrum can 
help reveal the physical conditions of the inner region around 
accreting compact objects. The X-ray spectrum of Cyg X-l 
while in the hard (i.e., low) state can be roughly described 
by a power-law with a photon-index V ~ 1.65, modified by 
an exponential cutoff with an e-folding energy Ei ~ 150 keV 
(Ebisawa et al. 1996; Gierliiiski et al. 1997, and references 
therein). A spectral shape of this form is naturally explained 
by Comptonization of low-energy seed photons by a semi- 
relativistic corona (Sunyaev & Triimper 1979; Titarchuk 
1994). Below about IkeV, there is evidence for a soft-excess, 

* Present address: Department of Physics and Astronomy, Uni- 
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usually interpreted as thermal radiation from a cold accre- 
tion disc having a temperature kT BB ~ 0.1 — 0.3 keV (Balu- 
rinska- Church et al. 1995; Ebisawa et al.1996). In addition, 
the spectrum is usually found to contain a weak iron line fea- 
ture at « 6.4 keV and a slight hardening above lOkeV, often 
interpreted as being due to Compton reflection (Barr, White 
&; Page 1985; Ebisawa et al. 1996, and references therein). 

The combination of a Comptonization continuum and 
reprocessing features has been interpreted as being due to 
an accretion disc corona (ADC). The geometric configura- 
tion of the corona and the cold disc, however, is still un- 
clear. In most prior work, the geometry has been assumed 
to be a cold accretion disc embedded between two hot coro- 
nae, with a slab geometry (Haardt, Maraschi Ghisellini 
1996; Haardt et al. 1993, and references therein). Recently, 
however, evidence has been presented showing that ADC 
models with a slab geometry suffer from several problems, 
making them less likely to be the appropriate models for ex- 
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plaining the high-energy radiation of BHCs. Gierlinski et al. 
(1997) analyzed simultaneous Ginga-OSSE data of Cyg X- 
1 and found that the strength of the reflection component 
is too small to be consistent with a slab geometry. Using a 
non-simultaneous broad-band spectrum of Cyg X-l, we pre- 
viously argued that our self-consistent ADC models with a 
slab geometry (in which the coronal temperature is not a free 
parameter but is determined by balancing Compton cool- 
ing with viscous dissipation) could not have a high enough 
temperature* to explain both the observed power-law and 
exponential cutoff (Dove, Wilms & Begelman 1997a; Dove 
et al. 1997b; Wilms et al. 1997). Moreover, the luminosity of 
the disc expected due to reprocessing is comparable to the 
coronal luminosity, and therefore the predicted soft-excess 
component is much stronger than that observed (Dove et al. 
1997b; Gierlinski et al. 1997). 

Most of the problems cited above are due to the cold 
disc of the slab ADC model having a covering fraction of 
unity (i.e., all downward directed coronal radiation is re- 
processed in the cold disc). ADCs having a geometry with 
a smaller covering fraction have weaker reprocessing fea- 
tures, are less efficiently Compton cooled, and thereby al- 
low higher coronal temperatures for a given r. An exam- 
ple is a hot coronal sphere that is surrounded by an exte- 
rior cold, optically thick, geometrically thin, accretion disc 
(called the “sphere+disc geometry” henceforth), similar to 
the two-temperature disc model of Shapiro, Light man & 
Eardley (1976) and to advection dominated models (Abra- 
mowicz et al. 1995; Chen 1995; Narayan & Yi 1994). Based 
on non-simultaneous BBXRT, Mir-TTM, Mir-HEXE, and 
OSSE data, we have shown that the sphere+disc ADC model 
can explain the observed spectrum of Cyg X-l for energies 
in the ~ IkeV— lMeV range (Dove et al. 1997b), 

In this paper, we analyze the broad-band spectrum of 
Cyg X-l from 3.0 to 200 keV, using data from the Rossi 
X-ray Timing Explorer (RXTE). The simultaneity of this 
broad-band observation allows us to place stronger con- 
straints on the coronal parameters than was previously pos- 
sible. In §2, we describe the observations performed and the 
procedures used in the data analysis. In §3, we give the re- 
sults from spectral modeling of the Cyg X-l spectrum us- 
ing both “standard” simplified models as well as our self- 
consistent ADC models. In §4 the results of these fits are 
discussed. 


2 OBSERVATIONS AND DATA ANALYSIS 

RXTE observed Cyg X-l for a total of 22.5 ksec on 1996 Oc- 
tober 23 and 24. For the analysis presented here, we use data 
from the proportional counter array (PCA) and the high en- 
ergy X-ray timing experiment (HEXTE). The PCA consists 
of five Xe proportional counters with a total effective area of 
about 6500cm 2 (Zhang et ad. 1993), while HEXTE consists 
of two clusters of NaI(Tl)/CsI scintillation counters with 
a total effective area of about 1400cm 2 (Rothschild et et. 
1997). For both instruments, only data where the source 

* Formally, a high temperature is possible for a slab corona if 
the coronal optical depth is very low. However, such low optical 
depth coronae also produce spectra that are too soft to explain 
the observations of Cyg X-l. 
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Figure 1. Top: Comparison of the Crab Pulsar spectrum, taken 
by PCA, to the best fit power law model. Bottom: Ratio of the 
residuals of the model to the data (data/model). Figure by K. 
Jahoda, private communication. 


was observed at elevations higher than 10° above the space- 
craft horizon were used. Due to missing standard- mode PCA 
data, we have used only the first 10 ksec of the observation. 
This duration is sufficient for the spectral analysis, as the 
PCA count-rate is high (« 4300 cps, with a background of 
» 150 cps). 

For the extraction of the PCA data, we used the stan- 
dard RXTE ftools, version 4.0. We used version 1.5 of the 
background estimator program, in which activation due to 
the South Atlantic Anomaly and an estimate for the internal 
background based on the measurement of the rate of Very 
Large Events within the PCA are taken into account. To this 
estimation, the measured diffuse X-ray background is added 
by the program. The comparison of the modeled background 
spectrum to the background measured by the PCA during 
the Earth occultation phases of the observation (defined by 
an elevation < —5°) indicates that the background model 
count-rate has an uncertainty of 5 to 10% (i.e., 0.3% of 

the total count rate). The use of the background model is 
still preferable to using the Earth occultation data because 
the latter method cannot estimate the variation of the back- 
ground due to the variation of the rigidity of the Earth’s 
magnetic field during on-source time intervals. 

For the spectral analysis of the PCA we used ver- 
sion 2.2.1 of the PCA response matrix (Jahoda, priv. comm.; 
for an overview of the PCA calibration issues see Jahoda et 
al., 1996). The PCA response matrix is now fairly well un- 
derstood. An analysis of the residuals of a standard power- 
law fit to the spectrum of the Crab pulsar indicates that the 
uncertainty of the calibration in the vicinity of the Xenon 
L edge (around 5.5 keV) is about 2%. This makes a precise 
determination of the parameters of the Fe Ka line difficult. 
In addition, the Crab fit indicates a calibration uncertainty 
of about 5% above « 30keV, as shown in Figure 1. Since 
HEXTE provides reliable data for energies £ 20keV we ig- 
nored the PCA data above 30keV. To be conservative, we 
ignored the first four PCA channels (i.e., we only included 
PCA data above 3keV). To account for the additional un- 
certainties in the PCA response matrix, a 2% systematic 
error was added to all PCA data. 

Software provided by the HEXTE instrument group was 

1-1 I ■ 
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HEXTE Observation of the Crab Nebulo/Pulsor 
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Figure 2. Top: Comparison of the Crab Pulsar spectrum, taken 
simultaneously by HEXTE Clusters A and B, to the best fit bro- 
ken power law model (Ti = 2.045 ± 0.002, T 2 = 2.140 ± 0.009, 
and the break energy E% = 57.2 ± 3.0keV). Bottom: Ratio of the 
residuals of the model to the data (data/model) (Rothschild et 
al. 1997, Fig. 11). 

used for the extraction of the HEXTE data and subsequent 
dead- time corrections. Since HEXTE is source-background 
swapping, the background is not a problem in the data anal- 
ysis. We used the HEXTE response matrices released 1997 
March 20. Only data above 20keV were used due to the 
uncertainty of the response matrix below these energies. At 
high energies, the spectrum was cut at 200 keV. To ensure 
good statistical accuracy above 50keV, the spectrum was re- 
binned by a factor of 3 for channels between 50 and 100 keV, 
and by a factor of 10 for higher channels. The spectrum of 
the Crab Pulsar, taken by HEXTE, is shown in Figure 2. 

3 SPECTRAL ANALYSIS 
3.1 Standard Models 

Spectral fitting was first performed using the following stan- 
dard models: a power-law, a power-law with an exponential 
cutoff, a power-law with an exponential cutoff plus a cold re- 
flection component, and thermal Comptonization models. In 
addition, we added a Gaussian line (with energy and width 
fixed to 6.4 keV and 0.1 keV, respectively) to several of these 
models. We fixed the low-energy absorption to an equivalent 
cold Hydrogen column of N H = 6 x 10 21 cm -2 , the value sug- 
gested by the soft X-ray spectrum and interstellar reddening 
measurements of HDE 226868 (Bahicidska & Hasinger 1991; 
Wu et al. 1982). The results of the spectral fits are given in 
Table it. 

For energies E 10 ke V, a power-law with an exponen- 

tial cutoff provides a good description of the data, as also 
shown in Figure 3. The value for the e- folding energy, £f, 

t To account for the known discrepancies in the relative nor- 
malization of PCA and HEXTE, we introduced a multiplicative 
constant in the spectral models that represents the relative nor- 
malization between the data-sets. The relative normalization be- 
tween HEXTE and PCA was always found to be 0.74 ± 0.01 and 
is therefore not listed in Table 1. 
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Figure 3. Comparison of the RXTE spectrum of Cyg X-l to the 
best-fit power-law with an exponential cutoff model. This model 
was fit to only the 10 - 200 keV data, yielding = 171/149, 
The photon power-law index T = 1.46 and the e-folding energy 
Ei = 168 keV. A soft-excess below lOkeV is clearly evident. 


is well constrained due to the quality of the HEXTE data. 
The analysis of the residuals of this fit indicates the pres- 
ence of a soft excess at energies 8keV, as shown in Fig- 
ure 3. Adding a black-body component with a temperature 
kT h B as 1 keV and a Gaussian component, with an equivalent 
width EW= 46.4 eV, to the model improves the quality of 
the fit, resulting in a good description of the data over the 
whole energy range (x?ed = 170/166), as shown in Figure 4. 
We note that this black-body component should not be in- 
terpreted literally as being due to disc emission; it is simply 
added phenomenologically to measure the magnitude of the 
data’s deviation from a power law. We could have instead 
used a disc blackbody, Gaussian, or bremsstrahlung model 
with roughly equal success. 

The addition of a reflection component to this model 
does not improve the quality of the fit, as the best-fit value 
for the fraction of incident radiation that is Compton re- 
flected by cold matter is / ^ 0.02. Even though we only 
considered PCA data in the 3.0 keV — 30keV range, the 
power-law index of the continuum is well constrained by 
the HEXTE data, and there is no evidence that the spectral 
slope of the data for the two instruments disagree within the 
energy range of the overlap. The PCA residuals between 20 
and 30 keV seen in Figures 3-6 are consistent with those 
seen for power-law fits to the Crab, as shown in Figure 1. 

We believe that the measured spectrum parameters, 
which represent a harder power law and a weaker reflection 
component than previously found by others, is mostly due 
to our fitting the entire 3 — 200 keV spectral band. Previous 
observations using different instruments have had relatively 
poor spectral coverage in the energy range ^ 30 — 100 keV. 
We do not believe that our results are heavily influenced by 
any remaining uncertainties in the PCA response matrix. 
To support this claim, if we restrict the data to 3 — 30 keV 
and 100 — 200 keV (roughly covering the energy range of 
the Ginga + OSSE observations) we do find best-fit param- 
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Gaussian + Blackbody + Cutoff Powertaw 



Figure 4. Comparison of the RXTE spectrum of Cyg X-l to 
the best-fit power-law with an exponential cutoff model. The 
photon power-law index T = 1.45 and the e-folding energy is 
E{ — 164 keV. A blackbody component, with a temperature 
kTsB = 1.1 keV and a relative flux of £ 5% } and a Gaussian 
component, with an energy fixed at 6.4 keV and an equivalent 
width EW= 46.4 eV, were added to fit the soft-excess. For this 
fit, = 170/166. 


Reflection Model 



Figure 5. Top: Comparison of the reflection model (Magdziarz 
& Zdziarski 1995) to the PCA data alone. Here, x? e d = 143/73. 
Adding a Gaussian component, with a line energy Et = 6.4 keV, 
a line width AE = 64 eV, and an equivalent width EW r = llOeV, 
results in a slightly better fit (Xr ed — 107/72), and the other best- 
fit parameters do not change appreciably. Bottom: Ratio between 
data and model. 

eter values similar to those previously found for Cyg X-l. 
Specifically, for this restricted data range, we find that a re- 
flected power law with F = 1.81 ±0.01 and covering fraction 
/ = 0.35 ±0.02 yields x?ed - 157/90 = 1.75. Adding back 
the HEXTE data in the 30— 100 keV range, significantly re- 
duces the quality of the fit (x?©d = 8.7, 77 dof), as shown in 
Figure 5. Being able to fit a reflected power-law with a large 


disc covering fraction is thus seen to be partly a consequence 
of ignoring the data in the ~ 30 — 100 keV range. 

It is this 30 — 100 keV energy range that is crucial for 
constraining a power-law index since these energies are un- 
contaminated by reprocessing or reflection features, i.e. , the 
fined power law in this energy range is identical to the intrin- 
sic power law prior to reflection We therefore postulate that 
the observed 10-30 keV “hardening” feature typically found 
by others appears to be the beginning of a hard power-law 
that continues to ~ 100 keV, and that the ~ 1 — 10 keV por- 
tion of the spectrum is softer due to contamination from the 
thermal excess. This contamination could be due to Comp- 
tonization of the thermal radiation emitted by the cold disc 
with m effective blackbody temperature fcT B B ~ 150 eV. 

We emphasize that this result does not indicate that 
there are no reflection features in the observed spectrum 
of Cyg X-l. As we discuss below, our ADC model with a 
sphere±disc geometry is able to describe the observation, 
and this model does include reflection. The reflection fea- 
ture predicted from our model, due to reprocessing from 
the Comptonized continuum (which is not a pure power-law 
due to the contribution from the seed photons, i.e., ther- 
mal emission of the disc), differs from the feature predicted 
by phenomenological reflection models such as the popular 
pexrav model (Magdziarz & Zdziarski 1995), in which a pure 
power law is incident onto a uniform slab. 

Although the residuals of our power-law fits clearly in- 
dicate the presence of an Fe Ka line, the ~ 2% uncertainty 
in the PCA response matrix at 6 keV make a determination 
of the line equivalent width problematic. The division of the 
Cyg X-l spectrum by a spectrum of the Crab pulsar (which 
is assumed to be a pure power-law) also shows that the Iron- 
line feature is present in our data. Due to different detector 
gains between the two observations, however, we cannot use 
the divided spectrum to estimate the equivalent width of 
the line. We did include a Gaussian line in most of our fits 
presented in Table 1. We fixed the line energy at 6.4 keV 
and the line- width to 100 eV [which was the typical value 
for the ASCA observations (Ebisawa et al. 1996); For all 
models, thawing the line- width resulted in a best-fit width 
of ~ 600 eV.J. Fitting only the line amplitude, the fits typ- 
ically yielded an EW a 60 ± 35 eV. Although undoubtedly 
influenced by the uncertainties in the response matrix, this 
EW does indicate a contribution from an intrinsic line and 
is probably in agreement with the results based on four days 
of ASCA observations presented by Ebisawa et al. (1996), 
who found a comparably weak Fe line with an equivalent 
width EW^ 30 eV. 

Thermal Comptonization is a more physical model for 
the observed energy spectrum of Cyg X-l. Applying the 
model of Titarchuk (1994) to the data, we find that a ther- 
mal Comptonization spectrum, resulting from a spherical 
geometry with r « 3.6 and kT c » 40keV, with an addi- 
tional soft component, can give an acceptable description of 
the data, although the model underestimates the flux at en- 
ergies 150 keV. The semi-relativistic, optically thin model 
(Titarchuk 1994), for either a slab geometry or a spherical 
geometry was unable to explain the data. 

Although the traditional models can be reasonably suc- 
cessful in describing the observed broad band spectrum of 
Cyg X-l, they do not yield a physical interpretation of the 
emitting mechanisms responsible for the production of the 

n! I 
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high-energy radiation. There need not be a one-to-one cor- 
respondence between a phenomenological model component 
and a physical interpretation. We therefore apply our self- 
consistent ADC models to the observed data. 

3.2 Accretion Disc Corona models 

For both a slab ADC model and the sphere-fdisc ADC 
model, we have computed grids of spectra using a non- 
linear Monte Carlo scheme based on the code of Stem et 
al. (1995) and modified by Dove et al. (1997a, b). The free 
parameters of the model are the seed optical depth r e , (the 
optical depth of the corona excluding the contribution from 
electron-positron pairs) and the heating rate (i.e., the com- 
pactness parameter) of the ADC. For a given compactness 
parameter, the temperature structure of the corona is deter- 
mined by balancing Compton cooling with heating, where 
the heating rate is assumed to be uniformly distributed. The 
e~e + -pair opacity is given by balancing photon-photon pair 
production with annihilation. Reprocessing of coronal radi- 
ation in the cold accretion disc is also treated numerically 
(Dove et al. 1997b)^. The outer radius of the cold disc is 
assumed to be five times the radius of the coronal sphere 
(a r = Rd/Rc = 5), although the results are very insensitive 
to the value of this ratio for a r £ 3 (Dove et al. 1997b). 
The model spectra have been implemented into the data re- 
duction software XSPEC (version 10.0) (Amaud 1996) for 
use in spectral fitting. The grids of model spectra and the 
interpolation routines, based on Delaunay triangularization, 
are available upon request. 

Slab-corona models do not result in good fits, as the 
predicted spectra are always much softer than the observed 
spectrum (x?ed ~ 80). This result is consistent with our pre- 
vious findings based on non- simultaneous data (Dove et al. 
1997b). The reason that the slab ADC models always predict 
a spectrum softer than observed is that there is a maximum 
coronal temperature for a corresponding total coronal opac- 
ity. As discussed by Dove et al. (1997b), no self-consistent 
slab model can have both a high enough temperature and a 
high enough opacity to yield a Compton- y parameter large 
enough such that the Comptonized spectrum is as hard as 
that of Cyg X-l. 

The sphere-fdisc ADC model does provide a good de- 
scription of the data (cf. Table 1 and figure 6)§. The formal 
X 2 -value of our best-fit model (*? e d = 155) is larger than 
the values found for some of the phenomenological models 

J One improvement to the sphere -f- disc model, not discussed in 
Dove et al. (1997b), is the treatment of thermal radiation emitted 
by the accretion disc. Currently, the temperature of the accretion 
disc is assumed to be Tub oc (R/Rc)~ 3 l 4 , where R is the disc 
radius and Rc is the radius of the coronal sphere. The implicit 
flux of the disc (i.e., the radiation emitted due to viscous energy 
dissipation and not due to reprocessing of coronal radiation) is 
given by F(R) = cT 4 (R ), and the spectral shape is determined by 
the local disc temperature (a superposition of thermal Planckian 
distributions). Reprocessing is treated locally and is added to the 
intrinsic thermal flux from the disc to yield the total flux. 

§ Note that, even though the free parameters of each grid of mod- 
els are the coronal seed opacity and the coronal compactness pa- 
rameter, for convenience we list the corresponding overage coronal 
temperature and total optical depth. 
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Figure 6. Top: Best-fit sphere-j-disc fit to the PCA and HEXTE 
data. Bottom: Ratio between data and model for the best-fit 
sphere-hdisc model. Deviations less than 2% are within the range 
of uncertainty of the PCA response matrix. Deviations in the 
20 — 30keV range are similar to those seen with best-fit models 
of Crab nebula observations (Jahoda, private communication). 


discussed in §3.1. Considering that our model has only three 
parameters (seed optical depth, coronal compactness, over- 
all normalization, plus the “hidden” parameter 7 b b, which 
is fixed at 150 eV during the fitting procedure), and consid- 
ering that our model is physically self-consistent, the level 
of agreement with the data is quite good. The major dis- 
agreement between the model and the data occurs in the 
« 5 - 10 keV range, where the residuals are as large as 7%. 
Below we discuss several ways in which slight modifications 
to our model might result in improved fits to the data. 

Contrary to the results of Gierlinski et al. (1997), we 
were able to fit the high-energy tail with a single Comptoniz- 
ing component. This result might be due to the ~ 15% vari- 
ation of the coronal temperature (due to the non-uniform 
radiation field and the corresponding non-uniform Comp- 
ton cooling rate) in our models. We also tried sphere-fdisc 
models in which the inner temperature of the disc is 300 eV 
or 800 eV, the latter temperature roughly being consistent 
with the best-fit black-body temperature found in the previ- 
ous section. However, since sphere-fdisc models calculate the 
flux emitted by the disc self-consistently, both of these mod- 
els predict much more of a soft-excess than observed. This 
result reinforces our conclusion that the IkeV black-body 
component that was added to the exponentially truncated 
power-law model should not be interpreted literally as disc 
emission. 
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4 DISCUSSION AND CONCLUSIONS 

We have applied a variety of spectral models to an RXTE ob- 
servation of Cyg X-l. We find that the observed spectrum in 
the 3.0— 200 keV range is well described by a power-law with 
a photon index T = 1.451® ^ modified by an exponential 
tail with an e-folding energy Ef « 160 keV, and a soft excess 
that was modeled as a black-body (although other broad dis- 
tribution models are capable of explaining the soft-excess) 
having a temperature fcT BB ~ 1.2 keV. The measured value 
of T is lower than those found in previous broad-band analy- 
ses, which find aT ~ 1.65 (Gierliriski et al. 1997; Dobereiner 
et al. 1994). In addition, the observed strength of the Comp- 
ton reflection feature is weak, as the best-fit covering frac- 
tion, determined by fitting the data with a power-law re- 
flection model (Magdziarz k Zdziarski 1995), is found to be 
a 0.2 when no soft excess component is included, and is 
found to be 0.02 when one is included. The latter value is 
considerably lower than the / ss 0.3 found by Gierlihski et 
al. (1997). 

Two weeks prior to our observation, Cygnus X-l tran- 
sited from a soft (‘high’) state to a hard (‘low’) state. For a 
description of the properties the Cyg X-l soft state, see Cui 
et ad. (1997). The proximity of our observation to this tran- 
sition might have influenced our results, leading to a harder 
spectrum and weaker reflection components than what has 
usually been observed from this source. However, we believe 
that our best-fit parameters are also due to our fitting the 
entire 3 — 200 keV spectral band, as discussed in section 3.1. 

ADC models having a slab geometry are unable to ex- 
plain the observed hard power-law and the small amount 
of reprocessing, consistent with the results of Dove et al. 
(1997b), Gierlinski et al. (1997), and Poutanen, Krolik k 
Ryde (1997). In contrast, the sphere+disc ADC models 
provide a good explanation of the data, a result in agree- 
ment with (Gierlinski et ad. 1997). Furthermore these models 
are physically self- consistent Our sphere+disc models seem 
to under-produce “reprocessing features,” especially the Fe 
line, even though reprocessing is self- consistently included 
within these models. For this paper, the ratio of the accre- 
tion disc radius to the coronal sphere radius, a was set to 
five. As shown by Dove et al. (1997b), for <3^5, the geo- 
metrical covering fraction of the accretion disc, as seen from 
the surface of the corona, is /g « 0.30. If these residuals, 
which represent less them 1% of the total flux from Cyg X-l, 
are due to underestimating reprocessing, several possibili- 
ties come to mind. First, our model has a sharp transition 
between the cold (flat) disc and the corona. It is possible 
that these two regions overlap to some extent, and thereby 
produce stronger reprocessing features. Poutanen, Krolik k 
Ryde (1997) suggest that these regions overlap to a large 
extent during the high state, which ended only two weeks 
prior to our observation. Another possibility is that we are 
observing a hot transition layer between the disc and corona. 
Our disc also remains flat out to its outer edge. If the disc 
flares, which is likely due to X-ray heating, then reprocessing 
features will also be enhanced. Finally, the iron abundance 
in the disc of Cyg X-l could be higher than the solar value 
that we used in our model. 

Due to the data’s broad energy range, the physical prop- 
erties of the corona are well constrained. However, it is pre- 
mature to claim that the sphere+disc configuration is indeed 


the appropriate geometry. Other, albeit unknown, geome- 
tries in which a small fraction of coronal radiation is repro- 
cessed by the disc, may also be able to describe the data. 
The spectral features due to the reprocessing of coronal ra- 
diation in the cold disc (i.e., the Fe Ka fluorescence line, 
the Compton reflection “bump,” and the soft-excess due to 
thermalization of the coronal radiation) occur for energies 
in the ~ 0.1 — 20keV range, and different geometries will 
predict slightly different reprocessing features. Since we had 
to ignore data below 3keV, we were not able to constrain 
the model parameters which deal with the soft-excess due 
to thermal radiation emitted by the cold disc, nor can the 
strength or width of the iron line be directly measured (due 
to the uncertainty of the PCA response matrix at 5.5 keV). 
Therefore, more low energy data are needed to further con- 
strain ADC models for Cyg X-l. 

Additional constraints to the ADC models can come 
from measurements of time-lags and the temporal coher- 
ence between several energy bands from the source. These 
measurements can be compared with the time-lags and co- 
herence function predicted using the physical parameters of 
the geometry found from spectral fitting (Vaughan k Nowak 
1997; Nowak k Vaughan 1996). Only a geometry in which 
both the spectral and the temporal data can be explained 
should be considered a valid candidate for Cyg X-l. We are 
planning to use such an approach in a forthcoming paper. 
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Table 1. Results of spectral fitting to Cyg X-l. 


Model 
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A PL 

E, 

kT B b 

A bb 

/ 

kT c 

T? 

Ac 

Al 

X 2 /dof 




[keV] 

[keV] 

10- 2 


keV 



10- 3 


Pi 

1.70 

1.78 


. , . 

. . . 



, , . 

. . . 


2315/171 

plexp 

1 A3+° 01 
AO,5 -0.01 

1 40+O.O3 
a * 4C -0.04 

364 1 3 ? 

. . . 


. . . 


. . . 

. . . 


806/169 

plexp+bb 

1 4C+0.01 
1,4 -0.02 

0 91+ 0 03 
u * yi -0.04 

1621? 

1 2+0.0 
A -0.1 

9 9 + ° 1 
*“*—0.2 

. . . 



. . . 


173/167 

plexp+bb+g 

1 45+° 02 
A,4D -0.0l 

0 91 +0.04 
u ' yi — 0.01 

1641? 

1 1 40 ' 1 
A - A -o.o 

2 l +0 2 
A — 0.1 



. . . 


3-3liI 

170/166 

pexrav 

1 71 + 002 
* a - 0.02 

1 67^"° 05 
A - D '-o.oe 

924l«° 0 



0 20^°* 05 
u * iju -0.05 





752/168 

pexrav +bb 

1 45+ 0 02 
A ' 4D _0.01 

n qo+o Q6 

164li° 0 

1 2 +0 1 
1 -0.1 

9 9 + ° 1 
* ‘*—0.2 

0 oi +00 ° 

V * U1 -0.0I 





173/166 

pexrav -{-bb-l-g 

1 44+ 0 02 

ri oq+0.04 

9-Oy_0.03 

16ll? 

1 2 +01 
A **-01 

1 9+° 1 
1,y -0.2 

n nn+o 02 




39l|;| 

156/165 

comptt+bb 




1 l +01 
A,1 -0.1 

lA -o'\ 

. . . 

40lJ 

3-61S.1 

O47+OO 1 

U ’ q/ -0.02 


244/167 

comptt+bb+g 




1 0 +0 1 
A ‘ -0.1 

1-aiSi 


40+ * 

36 ±ol 

0-47l®« 

™±ll 

232/166 

s+d 

... 

... 




... 

87l? 

i-6iS:i 

8.07l®“ 
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Nh was fixed at 6 x 10 21 cm“ 2 . Uncertainties given are at the 90% level for one interesting parameter (Ax 2 = 2.7). pi: power-law 
with photon-index T and normalization Apl (photon-flux at IkeV); plexp: power-law with exponential cutoff [Apl E~~ r exp(— E/Ef)]\ 
bb: black body with temperature fcX B b; the normalization A bb is defined in units of L^/R^ 0 where L39 is the luminosity in units of 
10 39 erg/s and Rio is the distance in units of lOkpc; pexrav: power-law with exponentional cutoff reflected off cold matter (using Green’s 
functions of Magdziarz & Zdziarski 1995), / is the ratio between the incident and the reflected flux; comptt: Comptonization spectrum 
after Hua & Titarchuk (1995) and Titarchuk (1994), with a coronal temperature of kTc and the total optical depth seed-photon 
temperature is 10 eV; g: Gaussian line, with energy and width fixed to 6.4 keV, 0.1 keV, respectively. Al is the amplitude of the line 
in total number of photons/cm 2 /s in the line. A l = 3.3 x 10~ 3 corresponds to an equivalent width of & 45 eV. s+d: ADC model with 
sphere+disk geometry, total optical depth of the corona is Tp, the coronal temperature is kTc » the temperature profile of the cold disk 
is *T bb (R) = 150 (R/Rc)" 3 ^ 4 eV. 
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